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ABSTRACT 

Recent theoretical investigations have suggested that the formation of the very first 
stars, forming out of metal-free gas, was fundamentally different from the present- 
day case. The question then arises which effect was responsible for this transition in 
the star formation properties. In this paper, we study the effect of metallicity on the 
evolution of the gas in a collapsing dark matter mini-halo. We model such a system 
as an isolated 3a peak of mass 2 x IO^Mq that collapses at Zcoii — 30, using smoothed 
particle hydrodynamics. The gas has a supposed level of pre-enrichment of either 
Z = 10"'' ^0 or 10~^Zq. We assume that H2 has been radiatively destroyed by the 
presence of a soft UV background. Metals therefore provide the only viable cooling 
at temperatures below lO^K. We find that the evolution proceeds very differently for 
the two cases. The gas in the lower metallicity simulation fails to undergo continued 
collapse and fragmentation, whereas the gas in the higher metallicity case dissipatively 
settles into the centre of the dark matter halo. The central gas, characterized by 
densities rtn ^ 10'* cm~'^, and a temperature, T ~ 90 K, which closely follows that of 
the cosmic microwave background, is gravitationally unstable and undergoes vigorous 
fragmentation. We discuss the physical reason for the existence of a critical metallicity, 
Zcrit ~ 5 X IO^^Zq, and its possible dependence on redshift. Compared to the pure 
H/He case, the fragmentation of the Z = 10~^Zq gas leads to a larger relative number 
of low-mass clumps. 

Key words: cosmology: theory - early universe - galaxies: formation - stars: forma- 
tion - hydrodynamics. 



1 INTRODUCTION 

One of the grand challenges in modern astrophysics is to 
understand how the cosmic "dark ages" ended, and to eluci- 
date the nature of the first luminous objects (e.g., Rees 1999; 
Barkana & Loeb 2001). Recently, various authors have ad- 
dressed the properties of the very first generation of stars 
through numerical simulations of the collapse and fragmen- 
tation of primordial clouds, consisting of dark matter (DM) 
and metal-free, pure H/He gas (Abel et al. 1998; Bromm, 
Coppi, & Larson 1999, 2001; Nakamura & Umemura 1999; 
Abel, Bryan, & Norman 2000). In the context of hierarchical 
models of structure formation, these so-called Population HI 
stars are expected to form in halos of mass ~ 10^ - 10** Mq, 
collapsing at 2; ~ 20 — 30 (Tegmark et al. 1997) . The three- 
dimensional numerical investigations have suggested that 
Population III stars were rather massive, perhaps even very 
massive with M. > IOOMq (Bromm et al. 1999, 2001; Abel 
et al. 2000). Based on two-dimensional studies, Nakamura & 
Umemura (2001) have hypothesized a bimodal initial mass 
function (IMF) for the first stars, with one mode favoring 



very massive stars, and the other leading to 'normal' stars 
of mass ~ IMq. The existence of a high characteristic mass 
scale has a robust explanation in the microphysics of cool- 
ing due to H2 which is the main coolant in the absence of 
metals. 

Although important uncertainties remain, mostly de- 
riving from the poorly understood physics of the protostel- 
lar feedback, it seems plausible that the formation of Pop- 
ulation III stars was qualitatively different from present- 
day (Population I) star formation. The question, therefore, 
arises: How did the transition in the star formation proper- 
ties take place, and what is the physics responsible for it? 
In this paper, we investigate the arguably most important 
effect: the presence of a trace amount of heavy elements. 
These metals were likely injected into the gas during a pre- 
enrichment event, since self-enrichment is rendered highly 
implausible due to the blowaway of any remaining gas after 
the first (Pop. Ill) stars have exploded as SNe, given the 
relatively small binding energy of Population HI DM halos 
(Ciardi et al. 2000). 

We here consider the limiting case in which all the H2 
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has been radiatively destroyed by the presence of a soft 
UV background. H2 molecules are very vulnerable to pho- 
tons with energies below the Lyman limit. Even before the 
universe has been reionized, a pervasive, soft UV radiation 
field could therefore have been established (e.g., Ciardi et 
al. 2000). Metals consequently provide the only cooling at 
temperatures T ^ lO^K. 

Other important roles in modifying the physics of star 
formation might be attributable to magnetic fields, or to 
the onset of turbulence. Both these effects are believed to 
be unimportant in the formation of the very first stars (e.g., 
Loeb 1998), but it will be interesting to explore their role in 
forming the second generation of stars in future work. 

Previous work on this question has focused on the chem- 
ical and thermal history of the gas, employing a simple one- 
zone model for the dynamical evolution of the collapsing 
cloud (e.g.. Silk 1977; Yoshu & Sabano 1980; Lin & Murray 
1992; Nishi & Tashiro 2000; Omukai 2000). Following the 
collapse to the point where the gas becomes opaque, Omukai 
(2000) argues that the resulting mass of the protostellar core 
is independent of metallicity. The final stellar mass, how- 
ever, is expected to have no direct relation to the initial 
(core) mass, and is mainly determined by the continued ac- 
cretion onto the protostellar core. Our approach, based on 
full three-dimensional numerical simulations of a collapsing 
primordial cloud, allows for the creation of sink particles 
and is thus well suited to address the complex physics of 
accretion. 

The paper is organized as follows. In §2, we present 
our numerical methodology. The choice and justification for 
the initial conditions, and the results of the simulation are 
discussed in §3. We summarize our findings, and address 
their implications in the final section. 



2 NUMERICAL METHOD 

The evolution of the dark matter and gas components 
is calculated with a version of TREESPIf (Hernquist & 
Katz 1989), combining the smoothed particle hydrodynam- 
ics (SPH) method with a hierarchical (tree) gravity solver 
(see Bromm, Coppi, & Larson 2001 for further details). Here, 
we discuss the additions to the code which are necessary for 
the investigation of low metallicity gas. These are a method 
to treat the radiative cooling of the gas, and a technique 
to create sink particles. The thermal evolution of the gas is 
governed by the equation: 

Dm _ P Dp r - A 

where D/Dt is the usual Lagrangian time derivative, P and 
p are the gas pressure and density, u is the specific internal 
energy (in erg g~^), and F and A are the contributions from 
radiative heating and cooling, respectively (in units of erg 
cm~^ s~^). The first term on the right-hand side describes 
adiabatic cooling due to expansion or heating due to com- 
pression. We now discuss the relevant radiative processes. 

We use the cooling function of Ricotti, Ferrara & 
Miniati (1997) which includes cooling by fine structure and 
metastable lines of C, N, O, Fe, S, and Si. Ionization equilib- 
rium is supposed to be maintained by cosmic rays that we 
assume to be associated with the heavy element production 
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Figure 1. Cooling function for gas which is enriched with a trace 
amount of metals. Shown is the cooUng rate (in units of erg cm^ 
s^^) vs. temperature, assuming a hydrogen density of Jijj = 1 
cm-3. Solid line: Z = IQ-'^Zq. Dashed line: Z = IQ-^Zq. 



by SNe. The (primary) cosmic ray ionization rate is scaled 
from the Galactic value, Ccfl ~ 1-8 x 10"^'' s~^ (McKee 
1995), by the factor Z/Zq, where Zq is the solar metallic- 
ity. The same scaling with Z is used to derive the dust-to-gas 
ratio, which we take locally to be equal to 1/160. Again, the 
rationale is that dust at high redshift can only be formed 
in the ejecta of primordial Type II SNe (Todini & Ferrara 
2001); hence, its production occurs simultaneously with the 
injection of heavy elements into the intergalactic medium 
(IGM). Dust, however, is only marginally important for the 
cooling (via electron recombinations on positively charged 
grains) at the low heavy element abundances we are inter- 
ested in. In Figure 1, we show the resulting cooling function 
for two different metallicities. 

Since radiative cooling to temperatures below that 
of the cosmic microwave background (CMB), Tomb = 
2.7 K(l + z), is thermodynamically not possible, we write 
for the cooling term 

A = A(T) - A(rcMB) . (2) 

For T < TcMB, radiative cooling consequently turns into 
heating. This approximate treatment ensures that T > 
TcMB, unless cooling proceeds via adiabatic expansion. At 
z ^ 30, gas temperatures are therefore limited to T ^ 90 K, 
justifying the neglect of cooling due to molecules, such as 
CO, which would become important only at lower tempera- 
tures. 

Due to the assumed complete destruction of H2 in the 
presence of a soft UV field, hydrogen molecules do not con- 
tribute to the cooling of the gas. Prior to reionization there 
are almost no ionizing photons, and we consequently ignore 
heating due to photo-ionization. 

We include heating due to the photo-electric effect on 
dust grains. This is consistent with our assumption that 
all the H2 has been destroyed by the presence of a soft 
UV background, since the photo-electric heating is mostly 
due to photons in the same spectral range. We write the 
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heating due to the photo-electric effect as F = an, with 
a — aQ{Z/ZQ). For the local value of the heating constant, 
we choose uq ~ 10~^'*e erg s~^, with a heating efficiency of 
e ~ 0.05 (Wolfire et al. 1995). As we discuss later, this effect 
is not very important for the evolution of the gas. 

We have devised an algorithm to merge SPH particles 
in high-density regions in order to overcome the otherwise 
prohibitive time-step limitation, as enforced by the Courant 
stability criterion. To follow the simulation for a few dy- 
namical times, we allow SPEI particles to merge into more 
massive ones, provided they exceed a predetermined den- 
sity threshold, nth — 10^ cm~^. More details of the merging 
algorithm are given in Bromm et al. (2001). 



3 THE SIMULATIONS 
3.1 Initial Conditions 

Within a hierarchical cosmogony, the very first stars are 
expected to form out of 3 — 4(t peaks in the random field 
of primordial density fluctuations. The early (linear) evolu- 
tion of such a peak, assumed to be an isolated and roughly 
spherical overdensity, can be described by the top-hat model 
(e.g., Padmanabhan 1993). We use the top-hat approxi- 
mation to specify the initial DM configuration, where we 
choose the background universe to be described by the 
critical Einstein-de Sitter (EdS) solution with density pa- 
rameters: Qdm ~ 0.95, Ob = 0.05, and Hubble constant 
h = _ffo/(100 km s~^ Mpc~^)=0.5. In this paper, we in- 
vestigate the fate of a 3a peak of total mass 2 x 10^ Mq, 
corresponding to 10^ M© in baryons which is close to the 
cosmological Jeans mass (Couchman & Rees 1986). On this 
mass scale, the standard Cold Dark Matter (CDM) scenario 
predicts a present-day r.m.s. overdensity of cro(M) ~ 16, 
with a normalization erg = 1 on the 8/i~^Mpc scale. Then, 
one can estimate the redshift of collapse (or virialization) 
from: 1 -I- Zcoii = 3cro(M)/1.69, leading to Zcoii — 30. 

Recent observations (see Barkana & Loeb 2001 and ref- 
erences therein) have presented evidence for a fiat universe 
with Qrn = 1 — = 0.3, thus rendering the EdS uni- 
verse increasingly implausible. At high redshift s {z ^ 10), 
however, the standard CDM model deviates only slightly 
from the more realistic ACDM case, and serves as a useful 
template for all hierarchical models of structure formation. 
The results presented in this paper mainly rely on the cool- 
ing physics of the gaseous component, and are not likely 
to change significantly, if these more realistic cosmological 
parameters are adopted. 

Our simulation is initialized at Zi = 100, by performing 
the following steps. The collisionless DM particles are placed 
on a cubical Cartesian grid, and are then perturbed accord- 
ing to a given power spectrum P{k) = Ak", by applying the 
Zeldovich approximation (Zeldovich 1970) which also allows 
to self-consistently assign initial velocities. The power-law 
index is set to n = — 3 which is the asymptotic small-scale 
behaviour in the standard CDM model (Peebles 1993). To 
fix the amplitude A, we specify the initial variance of the 
fiuctuations 

a^^AY^k" . (3) 
The summation is over all contributing modes, where the 



minimum wavenumber is given by the overall size of the 
Cartesian box, and kmax by the Nyquist frequency. Choosing 
erf ~ 0.1, the rms fluctuation at the moment of collapse is 

'^(^ = 30) = a, ^ 1 . (4) 

This choice ensures that the substructure develops on a 
similar timescale as the overall collapse of the background 
medium. Next, particles within a (proper) radius of Ri — 
150 pc are selected for the simulation. The resulting num- 
ber of DM particles is here A^'dm = 14123. Finally, the 
particles are set into rigid rotation and are endowed with 
a uniform Hubble expansion (see also Katz 1991). An- 
gular momentum is added by assuming a spin-parameter 
A = L\E\^^^/{GM^^^) = 0.05, where L, E, and M are the 
total angular momentum, energy, and mass, respectively. 

The coUisional SPH particles (A'^sph = 65536) are ran- 
domly placed to approximate a uniform initial density. The 
random sampling inevitably introduces shot-noise. The DM 
particles were set up on a regular grid speciflcally to avoid 
the unphysical shot-noise distribution which could mask the 
desired physical power spectrum. For the gas, however, the 
presence of the shot noise is not a big problem, since the 
gas mass is initially slightly smaller than the Jeans mass. 
Therefore, sound waves will efficiently wipe out all initial 
density disturbances. The SPH particles are endowed with 
the same Hubble expansion and rigid rotation as the DM 
ones. For the initial gas temperature, we adopt the value 
(see Tegmark et al. 1997): Tgas,i ^ 200 K. 

We carry out two simulations: Run A having a metal- 
licity oi Z = 10"^^© , and Run B one oi Z = W^Zq. 
These two values bracket those derived from QSO absorp- 
tion line studies of the low column density Lya forest at 
2: ~ 3 (Pettini 1999). They are also close to those expected 
to be associated with the reionization process. Otherwise, 
the two runs have identical initial conditions. 



3.2 Results 

The evolution of the collapsing clouds separates into two 
distinct stages. Leading up to the virialization of the DM, 
the two simulations behave very similarly, and the evolution 
is rather insensitive to metallicity. Once the gas has reached 
a roughly pressure-supported state in the potential of the 
DM halo, the subsequent evolution bifurcates, resulting in 
continued collapse and fragmentation for Run A, and in fail- 
ure to do so for Run B. We now discuss these two stages in 
turn. 



3.2.1 Virialization of Dark Matter Halo 

In response to the initially imprinted density fiuctuations, 
the DM develops a marked substructure. The gas does not 
'feel' the potential wells of these subcondensations, however, 
since it cannot cool sufficiently during these early evolution- 
ary phase. Instead, due to adiabatic compression, T oc n^^^, 
the gas reaches temperatures of ~ lO** K. At this point, very 
efficient cooling due to the excitation of hydrogen lines sets 
in (see Fig. 1), and maintains the gas at this temperature 
upon further compression. At the end of the virialization 
process, the gas has reached a state of rough pressure sup- 
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Figure 2. Run A: Morphology and gas properties at z = 30.5. 
Top row: The remaining gas in the diffuse phase. Left panels: 
Face-on view. Right panels: Edge-on view. The size of the box is 
30 pc. It can be seen that the gas has dissipatively settled into 
a central, disk-like configuration. Lower-left panel: Temperature 
vs. hydrogen number density. Lower-right panel: Jeans mass (in 
units of Mq) vs. hydrogen number density (in units of cm~^). 



port with typical gas temperatures close to the virial tem- 
perature 

GM rriH 



T 



5000 K 



(5) 



2A;b Rvir 

where Rvir — 100 pc is the virial radius, G is Newton's 
constant, fca Boltzmann's constant, and mn the mass of a 
hydrogen atom. To estimate the corresponding gas density, 
n„ir, we consider the baryonic Jeans mass 

and assume that it has to be approximately equal to the to- 
tal mass of the halo. Here, nn denotes the hydrogen number 
density. From imposing Mj ~ 2 x IQ^Mq, one then finds 
rivir — 10^'^ cm~'^. Heating due to the photo-electric effect 
is never important, neither during virialization, nor during 
the later, baryon-dominated stages. From here on, the fur- 
ther evolution strongly depends on the level of trace metal 
enrichment. 



3.2.2 Further Collapse and Fragmentation of the Gas 

We first discuss Run A. Cooling in this case is efficient 
enough to allow further collapse. In Figure 2, we show the 
gas morphology and the corresponding thermodynamic state 
at z — 30.5. As can be seen, the gas has dissipatively settled 
into a disk-like central configuration. This disk is horizon- 
tally supported by rotation. In the presence of a DM halo, 
one expects a contraction by a factor of 1/A ~ 20 for rota- 
tional support. The size of the disk, as shown in Figure 2, 
is in good agreement with this prediction. It is also evident 
that the disk is subject to a bar-mode (m = 2) instability. 
By examining the T — uu plane (lower-left panel in Fig. 2), 
one can see that the gas efficiently cools from ~ 10* K down 



Figure 3. Run A: Morphology at 2 = 28. Top row: The remaining 
gas in the diffuse phase. Bottom row: Distribution of clumps. Dot 
sizes are proportional to the mass of the clumps. Left panels: Face- 
on view. Right panels: Edge-on view. The size of the box is 30 pc. 
It can be seen that the gas has undergone vigorous fragmentation. 



to the value set by the CMB floor, T™™ ~ Tcmb = 86 K. 
During this rapid cooling, the gas remains rather smooth. 
The gravitationally unstable disk-material, however, subse- 
quently undergoes vigorous fragmentation (see Figure 3). 

We next turn to Run B, and to its rather different fate. 
In Figure 4, we show the situation at z ~ 28. In contrast to 
the corresponding situation for Run A (Figure 3) , no further 
collapse, and no fragmentation has occurred. Instead, the 
gas remains in pressure support, and lingers at temperatures 
close to Tvir. Only much later, as the result of slow but 
inexorable residual cooling, two high density clumps form 
at the very centre of the roughly spherical gas configuration 
(see Figure 5). Run B, therefore, presents an example of 
what Ciardi et al. (2000) have termed 'dark object', a halo 
which has failed to produce anything luminous in it. 

To summarize, in Run A, the gas becomes effectively 
self-gravitating, leading to continued collapse and fragmen- 
tation, whereas in Run B, the gas does not. To see this more 
clearly, we compare in Figure 6 the cooling timescale, tcooi — 
{nHkBT)/A, with the freefall time, tff ~ (GttihWh)"^'''^ 
(Rees & Ostriker 1977). Whether the gas can reach higher 
densities due to its self-gravity is decided at the end of the 
DM virialization process, with the gas having T ~ 5000 K 
and riH ~ 10^'^ cm~^. For these values, one has tcooi — tff 
for Run A, and tcooi > tff for Run B. Consequently, the gas 
in Run B is not able to cool any further and to reach higher 
densities. 

We now give an analytical estimate for the value of the 
critical metallicity, Zcrit, which discriminates between the 
two cases discussed above. Evaluating the cooling function 
(see Fig. 1) at T = 5000 K, we have 



A{Z) ~ 3 X 10"^°erg cm" 



The adiabatic heating. Tad 



10-4^(1 



P Dp 
p Dt 



can be written as 
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Figure 4. Run B: Morphology and thermodynamic properties of 
the gas at z = 28. The convention of Fig. 2 is adopted for the 
rows and columns. 



Figure 5. Run B: Morphology at 2 = 22.7. The convention of 
Fig. 3 is adopted for the rows and columns. 



^...10-ergc.-3s-(^)(Ji-,)^^ (8) 

Equating the two terms yields Zcrit — 5 x lO"*^©, in ac- 
cordance with the results of our simulations. 

We finally address the dependence of Zcrit on the red- 
shift of collapse, Zcoii, and the total mass of the halo, M. 
The gas temperature and density at the end of the DM viri- 
alization are expected to obey the relations T ~ T^ir oc 
M^/3(l -I- Zcoii) and nn tx pvir oc (1 + Zcou)^, respec- 
tively (e.g., Padmanabhan 1993). Assuming A (x T^''^ for 
10^ ^ T 10* K, inserting the virial relations into equa- 
tions (7) and (8), and normalizing to the case treated in this 
paper, one approximately finds 

This relation is only valid for halos with T^tr < 10* K. We 
plan to test this analytical prediction with future simula- 
tions. 



3.2.3 Comparison with the Z — Case 

Since Bromm et al. (1999,2001) have investigated the case 
of zero-metallicity gas with initial conditions that are other- 
wise very similar to those adopted here, it is instructive to 
compare our Z = 10~'^Zq simulation with the pure H/He 
case. In the previous papers, the formation of H2 was al- 
lowed for, whereas here we have considered the limiting case 
of maximum negative feedback which is supposed to destroy 
all H2. 

The first difference occurs during the virialization of 
the DM. In marked contrast to the simulations presented 
in this paper, the DM fluctuations do imprint their signa- 
ture onto the gaseous component if H2 is the main coolant. 
These DM induced condensations act as the seeds for the 
subsequent fragmentation of the gas. The basic reason for 
this difference is the boost in the formation of H2 due to the 
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Figure 6. The criterion for continued collapse. Freefall timescale 
{solid line) and cooling timescale {dotted symbols) vs. hydrogen 
number density (in cm"^). Timescales are in units of 10^ yr. 
Shown is the situation briefly after the virialization of the DM, at 
2 ~ 30.5. (a) Run A: The cloud succeeds in undergoing continued 
collapse due to efficient cooling, t^ool ^ 'tffy log "-H — 2.5. (b) 
Run B: Here, tcool > tfj sA, log nn — 2.5, and the gas fails to 
undergo continued collapse. 



density and temperature enhancement in the DM subcon- 
densations. Such a chemistry-related feedback is absent in 
the Z = 1Q~^Zq case. 

The second, crucial, difference derives from the fact that 
cooling remains efficient at densities nn ^ 10* cm"^, pro- 
vided the gas has been enriched to Z > Zcrit. As can be seen 
in Figure 2 (lower-left panel), the gas evolves isothermally 
up to the threshold density of nth = 10® cm"^. H2 cooling, 
on the other hand, saturates at nn ^ ricrit — — 10* 
cm""^, where Ucrit marks the transition from NLTE to LTE 
level populations. By comparing the distribution of Afj vs. 
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riH in Figure 2 (lower-right panel) to the equivalent Figure 
3 (panel (c)) in Bromm et al. (1999), it can be seen that 
there is no characteristic mass scale in the Z = 10"^^© 
case. Such a characteristic scale would correspond to a pile- 
up of SPH particles in the diagram due to long evolutionary 
timescales. In the pure H/He case, on the other hand, there 
is a preferred mass scale of M ~ IO^Mq. 

The overall range of clump masses, 10^ < Mci < 
is similar in the pure H/He and Z = 1Q~^Zq sim- 
ulations, although the relative number of low-mass clumps 
is significantly higher for the pre-enriched case. These low- 
mass clumps have masses close to the resolution limit of the 
simulation (Bate & Burkert 1997) 

Mres ^ Mb (^^^) - lOOM© , (10) 

where Nneigh — 32 denotes the number of SPH particles 
within a smoothing kernel. Obviously, clumps cannot have 
masses smaller than Mrea- The most massive clumps form 
with initial masses close to the Jeans mass, Mj ~ 75OM0, 
evaluated at the density, nn — 10* cm~'^, where the bulk 
of the freely-falling gas reaches the CMB temperature floor, 
Tmin ~ 86 K (see Fig. 2). 

Later on, clumps gain in mass through accretion from 
surrounding gas, and through occasional mergers with other 
clumps. An investigation of the possible subfragmentation of 
a clump lies beyond the resolution limit of our simulation. 
It is worth emphasising that the spectrum of clumps is dif- 
ferent from the stellar spectrum (i.e., the IMF), and it is 
not possible to directly derive the latter from the former. 
However, the qualitative trend that the relative number of 
low-mass clumps increases with metallicity is suggestive of 
the possible regulation of the IMF by the gas metal enrich- 
ment. To further constrain the clump mass spectrum, and 
eventually the stellar IMF, higher resolution studies will be 
necessary. 

The fact that the minimum temperature enforced by 
the CMB could imply larger Jeans masses at high redshift 
has already been pointed out by Larson (1998). The role of 
the CMB in setting the minimum temperature points to the 
possible importance of z^oii in setting an evolving mass scale 
for cosmic star formation (see Bromm & Clarke 2001 for a 
discussion of this effect). 



4 SUMMARY AND CONCLUSIONS 

We have discussed the evolution of primordial star-forming 
objects in the context of the standard CDM model for struc- 
ture formation. The studied clouds are pre-enriched to a 
level of 10~* and 1Q~^Zq, consistent with observations of 
the IGM, and collapse close to z^oii ~ 30. We find that 
the evolution proceeds very differently in the two cases. The 
gas in the higher metallicity run can efficiently cool, and 
consequently collapse to a disk-like configuration. The disk 
material is gravitationally unstable and fragments into a 
large number (~ 25) of high-density clumps. In marked con- 
trast to this case, the gas in the lower metallicity simulation 
cannot efficiently cool, and therefore remains in the post- 
virialization, pressure-supported state in a roughly spher- 
ical configuration. Such a system constitutes an example 
of a 'dark object', as discussed by Ciardi et al. (2000). 



These results indicate the existence of a critical metallic- 
ity, Zcrit ~ 5 X 1O~*Z0, below which the presence of heavy 
elements does not significantly influence the evolution of a 
primordial cloud. The value of the critical metallicity might 
be parameter dependent, being larger at lower z and in more 
massive halos. To account for the observed abundance pat- 
tern in metal-poor stars, Wasserburg & Qian (2000) have 
hypothesized an initial, prompt enrichment episode due to 
a generation of very massive stars. According to their model, 
these stars would have formed in gas with [Fe/H];$ —3. This 
nucleosynthetic evidence for a critical metallicity is consis- 
tent with our result which is based on the physics of star 
formation. 

Our investigation has treated the case of a low mass 
halo, of mass ~ 10^ Mq, and the question naturally arises of 
what happens in more massive ones. Is there still a critical 
metallicity, and if so, what would be its numerical value? 
In addressing these questions, one can argue that halos with 
virial temperatures in excess of~ 10* K are much less depen- 
dent on a level of pre-enrichment in order to form stars. The 
gas in these halos can initially cool through line-excitation 
of atomic hydrogen, and in the later stages of molecular 
hydrogen. Even if not present initially, or if completely de- 
stroyed in the shocks associated with the virialization pro- 
cess, H2 is expected to be efficiently formed in the post- 
shock regions of these systems (e.g., Shapiro & Kang 1987). 
Due to their larger binding energy, these halos can experi- 
ence multiple episodes of self-enrichment. Whether there is 
a strong change in the star formation properties, once this 
internal enrichment exceeds a certain value, cannot be an- 
swered without simulating the evolution of such a high-T^ir 
halo. 

Although there remain important uncertainties with re- 
gard to our findings, we conclude by discussing their po- 
tential implications. If indeed stars which formed from ex- 
tremely low metallicity gas, with Z < Zcrit, were predomi- 
nantly massive, this would entail short stellar-evolutionary 
timescales, tevoi ^ yr. Consequently, it is not expected 
to find any extremely low-metallicity stars which are still 
'alive' today. This prediction is consistent with observa- 
tions of metal-poor halo stars, having Z > 1Q~'^Zq (Beers 
2000). Assuming that massive Pop. Ill stars end their lives 
in massive black holes (MBHs), with typical mass ~ 10=* Mq, 
Madau & Rees (2001) have identified high-a peaks of mass 
~ lO^Af©, collapsing at z ~ 20 — 30, as the sites of their 
formation. As is suggested in this paper, MBHs could only 
have formed early in the history of the universe, out of gas 

with Z < Zcrit- 

If stars more massive than ~ 25OM0 buried their metal 
production in collapsing to a black hole (Fryer, Woosley, & 
Heger 2001), how did the IGM get enriched to a level of 
Z > Zcrit, enabling the subsequent formation of 'normal', 
low mass stars? The answer to this question has to await 
the precise determination of whether the majority of Pop- 
ulation III stars has masses below or above the threshold 
mass for the complete collapse into a black hole, and a bet- 
ter understanding of how the metals are expelled from the 
star forming halo (e.g., Madau, Ferrara, & Rees 2001). 
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